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In this work we study the Double Exchange-Holstein (DE-H) model with an electron-phonon inter- 
action 7 coupled to magnetism. The analysis is performed combining a mean-field approximation for 
the double exchange interaction and the Lang-Firsov transformation for the electron-phonon inter- 
action. Discontinuous magnetic transitions appear when the dependence of g with m is sufficiently 
large, resembling those experimentally observed in manganites. We observe that the characteristic 
resistivity peak that arises near the critical temperature appears for broad ranges of the system 
parameter values, unlike what occurs in a constant -7 model. 

PACS numbers: 71.38.-|-i , 75.10.-b , 75.30.Kz 



INTRODUCTION 

Perovskite manganites Lai-x-AxMnO^ {A = 
Ca,Sr,Ba) have been studied intensively since the 
discovery of the spectacularly large dependence of 
their resistivity with the applied magnetic field, a 
phenomenon called colossal magnetoresistance (CMR). 
This unusual high response to the application of a 
magnetic field occurs mainly for temperatures near the 
critical temperature of a ferromagnetic-metallic (FM) 
to paramagnetic- insulating transition (PI) In this 
range of temperatures, the curve of resistivity versus 
temperature exhibits a characteristic pronounced peak. 
It was found that the ferromagnetic-to-paramagnetic 
transition might be of first or second order [2] , depending 
on the composition of the material. 

The connection between ferromagnetism and metallic 
behavior can be explained with the Double Exchange 
(DE) model proposed by Zener Q. This model assumes a 
very large Hund coupling between electrons and localized 
spins, which causes a reduction of the effective electron 
hopping when the ion spins are disordered Q. Thus, 
the model predicts that the system is metallic in the fer- 
romagnetic state and an insulator in the paramagnetic 
one. However, the magnetoresistance effect that is ob- 
tained with this theory is very poor in comparison with 
the experimental observations, suggesting that the dou- 
ble exchange mechanism alone is not enough to explain 
the properties of the manganites Millis, Littlewood 
and Shraiman proposed that a minimal model for the 
manganites must include a term describing a coupling 
between electrons and lattice degrees of freedom. The 
effect of such interaction was considered by Roder, Zang 
and Bishop j[6[ , and by Millis, Shraiman and Mueller 
Q . In Ref. [0] the authors showed that the e-ph interac- 
tion gives a dependence of Tc with doping that is similar 
to the one observed experimentally. In Ref. f^, it was 



found that the e-ph coupling enhances the magnetoresis- 
tance effect and does provide a resistivity behavior that 
is in qualitative agreement with the experimental obser- 
vations. However, the resistivity peak near Tc is obtained 
with this theory only if the electron-phonon coupling 7 
satisfies the condition 1.08<7<1.2. A similar result 
was obtained by Verges et al using Monte Carlo simula- 
tions @]. The obtained phase diagram indicates that for 
the case n = 0.08, there is a FM to PI transition if 1.3 
< 7 < 1.6. 

In a recent paper, Dagotto et al studied the resistivity 
of the DE model including e-ph interactions using Monte 
Carlo computational techniques Q. In agreement with 
the previously mentioned works, it was found that the 
maximum of resistivity is observed in the DE-H model 
if the coupling 7 is fine tuned around some values that 
depend on electron density n. For instance, for n = 0.1, 
the maximum resistivity appears for values of 7 in the 
range (1.3, 1.6), in agreement with the results of [8|]. The 
inclusion of disorder in the DE-H model was considered 
by Kumar and Majumdar in [l3l and also in [§]. It was 
found that quenched disorder favors polaronic formation 
and the appearance of the resistivity peaks in the curve 
p{T). However, it is not clear whether the magnetoresis- 
tance effect induced by disorder reproduces the experi- 
mentally observed resistive behavior in all their aspects. 

In this work we propose an alternative modification of 
the model. We study an DE-H with an electron-phonon 
coupling 7 that depends on the magnetic ordering, in- 
stead of being constant as it is considered in the usual 
model. The variation of 7 with disorder has been a long- 
standing problem. Both theoretical and experimental re- 
sults show that 7 may have important changes induced 
by disorder 11, l^, 13] ■ It is therefore of great interest 
to investigate how the magnetic transitions are affected 
by the variations of the e-ph coupling. The work is orga- 
nized as follows. In section 2, after performing a Lang- 
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Firsov transformation of the Hamiltonian, we construct a 
free energy for the system following Kubo-Ohata theory. 
In section 3 we discuss the dependence of the electron- 
phonon coupling with the magnetization. In section 4 
we shall consider the effect of this dependence on the 
magnetic transitions and the resistivity. In section 5 we 
summarize our results. 



HAMILTONIAN MODEL 

The Hamiltonian of the Double Exchange-Holstein is 
given by 

H = Hz + Hh 

Hz = - ^ Uj {c\^Cja + h.c.) - ^ cTi • 

ij i 

Hh ^^lY.''U^-^b\ + h)+ojY.^b\h + \) (1) 

i i 

where is the electron hopping between the sites i and 
j, cl^ is the creation operator for the itinerant electrons, 
Jh is the Hund coupling, at is the spin of the itiner- 
ant electrons. Si is the spin of the Mn ions, 7 is the 
electron-phonon coupling and &J are the phonon creation 
operators. Here we consider 5* = 3/2, which is the to- 
tal spin of the three electrons of the Mn^^ ions. To 
study this Hamiltonian, we assume the usually consid- 
ered regime Jn/t — > 00. It has been shown by various 
authors (Q, 0)' that in this limit Hz reduces to 
the DE Hamiltonian: 



s 

z{X) = exp(A//5). (5) 

i=-s 

Since in the present model 7 depends on the magneti- 
zation, it must be also a function of the field parameter 
A. The value of try is also obtained averaging over the 
states of the dimer Substituting cTy by it averaged 
value cr, the Hamiltonian ^ becomes 

H= -criX) ^ ty(cLc,. + /i.e.)- 7(A) ^cLc..(6l-)- 6.) 

ij i 

+ u;J2{blh+^-). (6) 

i 

To treat the phonon dependent part of ^ we use the 
Lang-Firsov transformation H = e^He^^ , where S — 
ni{ai — a|) and g — '^/oj. With this transformation, 
the Hamiltonian takes the form 

H = ^a{X)e-^'J^uM.c,.e-^^-'''^e^^'-'"'^ + h.c.) 

i i i 

where the number operators for the electrons and 

ni — riif -|- rii^ . We approximate the wave function of the 
system as a tensorial product of waves function for the 
electrons and phonons 1-0) = 1-0)6 ® |V'o)p/i- Considering 
that \^o)ph is the vacuum state and averaging over this 
state, the Hamiltonian takes the form 



Hde = ^'^UjCTijicl^Cja + h.c.) 



(2) H= a(A)e-«^(^)^t.,(cU-. + /i.c.)-c..g2(A)^. 



where 



LJrp 



2S+1 



(3) 



Sji being the total spin of the subsystem formed by the 



ions at sites i,j and the electron, i. e 



•^0 



|Si-l-S,- -t-o-jl 



14 1 , , [31 ■ The value of ct^ must be obtained averaging 



over all the states of . To calculate different averages, 
we follow the mean-field approach of Kubo and Ohata 
and introduce an effective field A — h^ff/T that tends 
to order the ion spins 
obtained from 



1J|,[16|. The magnetization is 



2u^g\\)Y, 



(8) 



We note that the Lang-Firsov transformation introduces 
a Hubbard like attractive interaction which can promote 



However, we shall 



the formation of bipolarons [17 . 
consider the system is in the adiabatic limit lu ^ and 
thus we shall omit the last two terms of the above equa- 
tion. From ^ we can obtain a free energy F = E — TS 
for the system. The free energy of the electrons are given 
by 



f^e(A) 



„/3(M-efc(A))l 



2/3 



(9) 



m(A) = J2 il/S)exp{\l/S), 



i=-s 



(4) where Cfe = (T(A)e ^^'^^eofc, and eofe stands by the energy 
levels of the bare band. To obtain the free energy of the 
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whole system, we follow Kubo and Ohata lj| , adding to 
rie an entropy term corresponding to the spin ions. The 
free energy is thus given by 



T = n^iX) - 7Vr[lnz(A) - Am(A)]. 



(10) 



In order to determine the mean-field solutions at given 
temperature and doping level, we minimize J- with re- 
spect to the variational parameter A. We still need, how- 
ever, a dependence of g with A, or equivalently, a depen- 
dence of g with TO. 



COUPLING BETWEEN g AND THE 
MAGNETIZATION 



To justify the use of an electron-phonon parameter 
that varies with magnetization, we examine now the cal- 
culation of the electron-phonon coupling. With this pur- 
pose, we can expand the interaction between the electron 
and the ions Vd in the form 



where Jhiu ■ Si) is the spatial average of the Hund inter- 
action. The interaction electron-phonon takes the usual 
form, if we except the factor (1 — ^ J//((t • Si)). The usual 
electron-phonon interaction that does not depend on the 
magnetization is recovered when Jh = 0. In order to ob- 
tain the approximate dependence of g with to, we again 
follow the method used in 14 1 , Q and consider a dimer of 
two ions with one itinerant electron. If we suppose that 
initially the electron is at the ion i = 1 , in virtue of the 
strong Hund coupling the spin is aligned with the spin 
of this ion. When the electron jumps to the site i = 2, 
it will interact with the spin 82- In general di and S2 
will not be aligned. To obtain the value of Jh(cr • Si) 
we average Jh{cri • S2) over all the states of the total spin 
St — cr-l-Si-|-S2 in the presence of the field A. Due to the 
fact that a is aligned with Si, these two spins can be con- 
sidered as an unique spin with module S = S + ^. Since 
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S = 3/2 and a =1/2, {a- S2) =_(Si • S2)/4. To calculate 
this average we use the identity Si • S2 = 5(8^ 
Then, the average is explicitly obtained from 



-Si) 



= Yl K,:(r, -R°) QzV,yez(r, -R^ 



(11) 



where Qi measures the separation of the ion position Ri 
from it equilibrium value R^ i.e. Qi = Ri — R^. The first 
term in (fTTI) represents the potential that interacts with 
the electron in the absence of lattice deformation. Then, 
the information concerning the electron-phonon interac- 
tion is contained in the second term [l^ . We can separate 
the interaction electron-ion in the form T4i = Vsa + Vei , 
where Vsa corresponds to the Hund coupling and Vei in- 
cludes all other possible electron-ion interactions. Sup- 
posing that the position dependence of Vs^ is of the form 
Vsa = Vso-(ri — Ri)(— •/_ffcr • Si), the second term in the 
expansion pTjl may be put as 



Ve^2 Q»V,(T4i(i-, - R" ) - Jh<J ■ S.V,s(r, - R°)). 

ij 

(12) 

From the above expression, it becomes clear that Hund 
interaction can introduce large changes in the interac- 
tion electron phonon, principally when Jh is the strong 
coupling regime. To obtain a relationship involving to 
and g, we shall assume that the form of 14; (r^ — R?) 
and Vasi^j ~ R?) ^'^s similar. In the appendix, we 
show that this approximation is adequate under cer- 
tain conditions. Using this simplification we can write 



Vei(r, 



R" 



iVcrs{i^j — Ri'), where ^ is a constant. 



Inserting this equality in p2p . we obtain that 



(ai • S2) - 

E 51+52 sr^5T 
Sr = l/2 I^M=~St 



{SM\S^ 



s? 



St 



(14) 



where (5Af |S|, - Sf - S2)|S'Af) = {St{St + 1) - 82(82 + 
1) - 8i{8i + 1)), with 81 = 2 and ^2 = 3/2. Since 
TO and (tr • Si) are functions of A, we can obtain the 
relationship between these quantities. This dependence 
is shown in Fig. 1, and it can be accurately approximated 
as (a ■ Si) = ^8 m^-^. {a ■ Si) takes it maximum value of 
in the fully polarized state to = 1, when a and Si are 
aligned, and becomes zero in the disordered state to = 0. 
Then, from ([13]), the dependence of the electron-phonon 
coupling with to can be expressed as 



9 = 5o(l - am ■") 



(15) 



where a is a constant and go is the value of the e-ph cou- 
pling in absence of magnetization. From this expression, 
a can be expressed as a = {go — gm=i) / go, i-e. is the vari- 
ation of g between the paramagnetic state (to = 0) and 
the fully polarized state (to = 1) relative to go- We no- 
tice that in previous works Q, the notion of an effective 
e-ph coupling that changes with T was already employed. 
However, the variations of g considered here are not only 
effective but real ones. As we shall show later, the real 
variations of g introduces important modifications in the 
behavior of the system. 



NUMERICAL RESULTS 



Ve^2 - -(l-e^H(fT,.S,))^Q,V,Ks(rj-R°)), (13) 



The numerical solutions are obtained minimizing free 
energy subject to the conditions (|3|) and (fT5|) . In 
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Fig. 2 we show the magnetization as a function of T for 
a = 0.4, a doping level x = 1 — {rn) = 0.25 and different 
values of . Here a square DOS bare band of bandwidth 
D = 12t was assumed. In general, the effect of the e-ph 
coupling is to reduce the magnetic critical temperature, 
for either of the cases in which g is dependent or inde- 
pendent of 771. This occurs due to the narrowing band 
effect of the electron-phonon coupling, which reduces the 
efficiency of the double exchange mechanism. However, 
when g is m-dependent, the magnetic transition becomes 
discontinuous as long as a is above a critical value etc 
that depends on the value of go ( see Fig. 3 ). We notice 
that as a increases, both the critical temperature and the 
value of m may increase. This is because larger values of 
a reduces the value of g in the magnetic state, increasing 
the bandwidth and favoring the DE mechanism. On the 
other hand, when the value of a is increased, the critical 
value of (7o for the ocurrcncc of discontinuous transitions 
decreases. The minimal value of f/o for the occurrence 
of first order transitions goes to infinity as a — > 0, indi- 
cating that there are not discontinuous transition when 
a = Q { g independent of m) within our mean-field treat- 
ment. In Fig. 4 we show the phase diagram in the 
( 5o J Q! ) space. The labelled regions correspond to: 
I) continuous ferromagnetic-transitions, II) discontinu- 
ous magnetic-magnetic transitions and HI) discontinuous 
ferromagnetic-paramagnetic transition. The three cases 
are depicted in Fig. 3. Figs. 2, 3 show that the magnetic 
transitions obtained with an m-dependent g compares 
well with the one obtained in experiments when a > 
and go is strong. At this point we remark that the de- 
tailed form of the dependence of g on m does not affect 
significatively the system behavior. If we employ a sim- 
ple linear dependence of the form g = 170(1 — am), the 
phase diagram is nearly the same as those showed in Fig. 

4, although there are some changes in the critical tem- 
peratures. Although we have not calculated the value 
of a from first principles, it is worth mention that the 
first order transitions appear for any value of a > 0, if 
go is sufficiently large. Experimentally, it was found that 
the electron-phonon interaction may undergo variations 
of order 150 % [13] ■ If we assume this level of variations, 
the value of a may be taken between and 0.6. 

An important property of the manganites is the strong 
dependence of their resistivity with the application of a 
magnetic field near Tj.. In order to study if this effect 
appears in the present model, we examine the depen- 
dence of the bandwidth with an applied field. The av- 
erage physical magnetization in z direction is given by 
M = gs^BN{Sz), where N is the number of ions. Thus, 
the energy due to the magnetic interactions is given by 
Hm = Nhm, with h = gsHsSH. We then added the 
term Hm to the free energy (|18p and obtained the mean- 
field solutions. The variations of W are shown in figure 

5. As it can be seen, when the transitions are discon- 
tinuous (a = 0.4), there is a strong dependence of the 



bandwidth with H near Tc. Since W is proportional to 
conductivity, this means that there is a large increment 
of the conductivity with application of H. This effect is 
relatively weak when the transition is smooth (a = 0.1). 

In many works it has been reported that the first order 
transitions in manganites are accompanied by phase sep- 
aration. In order to analize this posibility we plotcd the 
free energy of the F and P solutions as a function of dop- 
ing for constant temperature. Tendency to phase separa- 
tion is estabhshed when the stability condition {§j^)t ^ 
is not satisfied. Each solution satisfies this condition 
separately, but at the point at which the two free energies 
cross, the first transition takes place and the condition 
is not satisfied. We then search for bi-phase solutions of 
the form: T = NiJ^i{Nei/Ni) + iV2J^2(A^e2/^2), where 
Ni are the number of sites ocuppied by phase i and N^i 
the numbers of electrons in the volume occupied by phase 
i. Thus, Ni+N2 = N and Nei+Ne2 = Ne, being Ne the 
total number of conduction electrons. jFi and J-2 are the 
free energies for each phase, showed in Fig. 6. Defining 
ai — Ni/N, zi — Nei/Ne and n — N^/N, we can express 
the free energy per site as: 

J^iT, n)/N = ziTi{T,nzi/ai) + il~zi)T2iT, itlill). 

1 — di 

(16) 

We then obtained the volume fraction occupied by each 
phase minimizing J- with respect to zi and ai , with the 
restriction that the densitiy of electrons rii — N^i/Ni 
must range between and 1 in each phase. The vol- 
ume fraction occupied by each phase is shown in figure 
7, from which we can see that phase separation takes 
place around the critical conditions for the first order 
transitions. 

As mentioned previously, the origin of the resistivity 
maximum that appears near the critical temperature is 
still not well understood. As already mentioned, this 
peak is obtained in the DE-H model [3]) 0,0, when the 
electron-coupling 7 takes values in a relatively narrow 
range. The calculations of the resistivity in these previ- 
ous works show that p decreases monotonically with T if 
7 < 1.2. On the other hand, if 7 > 1.6 , the resistivity 
rapidly increases with 1/T. This last behavior is prob- 
ably due to the formation of polarons, whose resistivity 
exhibit such a dependence with T. Although the spin 
scattering could contribute to develop a maximum at 
the resistivity of the polarons increase so rapidly that the 
peak is not observed. In order to estimate the variation 
of the resistivity due to polaron formation, we shall use 
the expression for ppoi in the strong coupling adiabatic 
limit {uj/t 0) 

Ppo,(r) = ATexp(-^) (17) 

where Ca ~ ep/2 is the activation energy f]\ and ep = g'^uj 
is the polaron binding energy. From the Kubo's formula 
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for resistivity it follows [20| that p ex Then, we can 

write A oc where A is the parameter appearing 

in (|17p . In various works it was found that this expres- 
sion agrees very well with the resistivity in manganites 
in the region T > [IH, [H,!!!!, revealing the pres- 
ence of small polarons in the paramagnetic phase. It is 
worth mention that in principle it is difficult to differ- 
enciate the polaronic resistivity from the resitivity of an 
ordinary semiconductor, which has a similar dependence 
to (fT7|) (without the prefactor T). However, in Ref. [2^ 
it has been shown that a band insulator theory cannot 
explain the differences between the gap value obtained 
from thermopower and the one obtained from resistivy 
measurements. On the other hand, these differences may 
be explained assuming polaronic interactions. 

We studied the variation of ppoi [T) using pT|) with the 
values of g and W obtained in the mean-field solutions. 
In order to consider a realistic value of w, we included the 
third term in ([8]) and minimized the free energy given by 



T = - A^T[ln z{\) - Ato(A)] - ug"^ (A) 



(18) 



In order to take realistic values of the parameters, we 
used the estimations of the activation energy and phonon 
frequency given in [i^l, Ca = 0.08 eV, w = 0.036 eV, 
which gives 50 = 2. On the other hand, the bare band- 
width has been estimated in Q to be of order 2.5 eV. 
Since in our model D = 12t, this gives cu ^ 0.2t. In Fig. 
8 the values of log{ppoi) for go = 2 a ~ 0.4 and uj = 0.2 
{t = l)are shown. From this figure we can see that a 
peak accompanied by a pronounced depletion of p is ob- 
served at T = Tc- This is caused by a rapid reduction of 
g ( in this case from g = 2 in the P state to g = 1.2 in 
the F state). As it can be seen from Fig. 8, relatively 
small changes in g produce drastic changes in ppoi. Mod- 
erate variations of 7 or other system parameters will not 
suppress the presence of these peaks. These appear in 
region III of Fig. 4. Of course, we have not considered 
the scattering of the electrons by the spins. We expect 
that this effect will superpose to the polaronic ones and 
contribute to increase the resistivity peak. 

We note that the reduction of t; allows a transition from 
small polarons (large coupling) in the P state to large po- 
larons (small coupling) in the F state. This transition has 
been reported by a number of groups Q , [2^ , [s^ . The 
variation of the resistivity for T < Tc also supports reduc- 
tion of the electron-phonon coupling in the ferromagnetic 
region. It has been shown that resistivity in the ferro- 
magnetic phase of the manganites can be well described 
with a dependence of the form p = po + AT^ + BT"^-^, 
where the first term is a constant, and the two last terms 
correspond to electron-electron interaction and electron- 
magnon scattering respectively [13] , 21 1 . A kind of tran- 



sition between two polaronic regimes, as the one found 
here, could explain the small influence of polaronic ef- 



fects in the resistivity below T^. We comment that we do 
not attempt to describe accurately the resistivity of the F 
phase with (fTT]) . because it is valid for the strong coupling 
regime. Thus, the resitivity of Fig. 8 describes the tran- 
sition between two phases with different e-ph coupling, 
both in the strong coupling case. It is worth mention 
that the form of this curve is similar to the one obtained 
in Loo.ssS'ro.isMnOa. At Tc — 240i^, this compound 
experiments a transition from a P insulator to an F in- 
sulator [iil (see inset in Fig. 40 of Ref. fl). 

A recent study of the quasiparticle excitation spec- 
trum of Lao.77Cao.23Mn03 using scanning tunnelling 
microscopy by S. Seiro et al revealed that the spectra 
presents a polaronic gap in the F and P phases. How- 
ever, the gap is reduced when the system enters in the F 
phase. Since the polaronic gap measured in the experi- 
ment is of magnitude Cp, the dependence of the gap with 
the temperature is compatible with the reduction of the 
e-ph interaction in the F state that is obtained in the 
present model. 



CONCLUSIONS 

In this work, we study a DE-H model with an electron- 
phonon interaction that depends on the magnetic order- 
ing. By a simple argument, examining the calculation 
of the electron-phonon coupling, we show that g is af- 
fected by the magnetic interactions. This occurs prin- 
cipally when the Hund coupling is strong, which is pre- 
cisely the regime that is considered to be relevant for 
the manganites. Introducing a field parameter to control 
the magnetization and making some assumptions about 
the spatial form of the electron-ion interaction, we cal- 
culated a dependence of g with the magnetization m. 
When a (~ dg/dm) is above a critical value that de- 
pends on (7o, the F-P transitions becomes discontinuous. 
This is not only because the DE mechanism and the po- 
laron formation are affected mutually in an indirect form, 
but also because there is a interplay between these two 
mechanisms. In the spin disordered state, the electron- 
phonon interaction is enhanced, favoring the appearance 
of small polarons. In turn, the polaronic effect reduces 
the electrons mobility decreasing the efficiency of the DE 
mechanism. This highly non linear effect induces a sharp 
magnetic transition. These first order transitions may 
occur directly between the F to P states or between two 
F states. The last case appears only in a narrow region of 
the parameter space and probably will not be present in a 
more refined treatment of the model than the mean-field 
approach considered here. On the other hand, we can 
expect that a crossover from smooth to sharp transitions 
will persists in region II if fluctuations are included. The 
abrupt transitions of region HI (Fig. 4) have a shape 
that resembles the magnetic transitions experimentally 
observed in the manganites. We estimated the changes 
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in the polaronic resistivity using the expression of ppoi 
in the adiabatic hmit. A strong peak in the resistivity of 
polaronic origin appears at the critical temperature when 
the transitions are abrupt. This is very robust and ap- 
pears in region III of the phase diagram of Fig. 4. This 
allows to explain the presence of the resistivity maximum 
in a series of compounds, for which 7 and other properties 
may be different. For some values of the parameters, the 
discontinuous transitions are accompanied by a crossover 
from large polarons in the FM phase to small polarons 
in the PI phase, in agreement with several experimental 
observations 23, 1^ 131 . In summary, the obtained re- 
sults show that the DE-H model with an electron-phonon 
coupling that varies with the magnetization reproduces 
some properties of the manganites that are still not well 
understood. It is not our intention to neglect the im- 
portance of other factors such as structural disorder or 
phase separation in order to explain the manganites be- 
havior, but instead of proposing another mechanism that 
could be important to complete the description of these 
systems. 
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Appendix 

The spin dependent part of the electron-ion interaction 
is an exchange interaction, given by 



J{R) 



Mr)c^2{r')V{r' - r)4>2{r')(l)i{r)drdr' . (19) 



The orbital functions 4'i{r) are centered on different 
points separated by a distance R. In order to study 
the dependence of this interaction with i?, we follow 
Campbell et al [3l| and assume that the e-e interac- 
tion is a screened potential of the form V{r' — r) — 



r' I ) , being r — r' the distance between the 



electrons. We also assume that the functions (jji are of 
the form 0(r) — Cexp(— rer). Here 1/a is the screening 
length and 1/k is the localization length. We evaluated 
J for different values of a, k and R. We obtained that 
when a and k are nearly equal, the variation of J{R) 
follows the variation of V{R). When a > k, J{R) fall 
more rapidly than V{R). In the case a < k, the opposite 
occurs. Then, the assumption that J{R) has the same 
spatial dependence than V{R), it is valid for some partic- 
ular physical conditions. On the other hand, it is worth 
mention that the overall results do not depend strongly 
on the functional dependence of g with m. Thus, this 
assumption was performed to allow an analytical calcu- 
lation, but does not introduce results which disappear if 
other dependence is adopted. 
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FIGURE CAPTIONS 

Figure 1: Dependence of {a ■ Si) with the magnetiza- 
tion. The continuous line represents the values obtained 
with (HI and (dH). The dashed line corresponds to the 
approximate expression (a ■ S^) — 3/4to^'^. 

Figure 2: Magnetization as a function of the temperature 
for go = 2 and a = 0.1,0.2,0.3,0.4. Temperature is 
measured in units of t. 

Figure 3: Magnetization as a function of the temperature 
for a — 0.4 and go = 1, 1.5, 2. 

Figure 4: Phase diagram of the DE-H model with g = 
.go(l — am^-^). The labelled regions represent: I) contin- 
uous magnetic transitions, 11) discontinuous F-F transi- 
tions and III) discontinuous F-P transitions. 



Figure 5: Percentual variation of the bandwidth with 
application of magnetic field, calculated as {W(H) — 
W{0))/W{H), for g = 2, a = 0,0.2 and h = 0.02. W{H) 
and M^(0) denote the bandwidths with and without mag- 
netic field. 

Figure 6: Free energies per site over the bandwidth, for 
c, = 2, a = 0.2 and ksT =0.05 . The solid line corre- 
sponds to the F state and the dashed line to the P one. 
For these parameter values, phase separation takes place 
for doping value ranges < a; < 0.3 and 0.7 < a: < 1. 

Figure 7: Volume fraction occupied by F (solid line) and 
P (dashed line) regions, ior g = 2, a — 0.2 and ksT 
=0.05 . 

Figure 8: Logarithm of the polaronic part of the resistiv- 
ity obtained from the expression ppoi{T) — exp(-^^), 
using the values of g and W obtained in the mean-field 
solutions for go — 2, oj — 0.2 and a = 0.4 (B is a con- 
stant). 
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